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Analysis of Some Singular Solutions 
in Fluid Dynamics* 

Zhouping Xint 

Abstract 

Studies on singular flows in which either the velocity fields or the vorticity 
fields change dramatically on small regions are of considerable interests in 
both the mathematical theory and applications. Important examples of such 
flows include supersonic shock waves, boundary layers, and motions of vortex 
sheets, whose studies pose many outstanding challenges in both theoretical 
and numerical analysis. The aim of this talk is to discuss some of the key issues 
in studying such flows and to present some recent progress. First we deal with 
a supersonic flow past a perturbed cone, and prove the global existence of a 
shock wave for the stationary supersonic gas flow past an infinite curved and 
symmetric cone. For a general perturbed cone, a local existence theory for 
both steady and unsteady is also established. We then present a result on 
global existence and uniqueness of weak solutions to the 2-D Prandtl's system 
for unsteady boundary layers. Finally, we will discuss some new results on the 
analysis of the vortex sheets motions which include the existence of 2-D vortex 
sheets with reflection symmetry; and no energy concentration for steady 3-D 
axisymmetric vortex sheets. 
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1. Introduction 



Many physically interesting phenomena involve evolutions of singular flows 
whose velocity fields or vorticity fields change dramatically. Shock waves, vortex 
sheets, and boundary layer are some of the well-known examples of such singular 
flows which provide better approximations for significant parts of the flow fields 
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near the physical boundaries, in the mixing layers and trailing wakes, etc., at high 
Reynolds numbers. The better understanding of the dynamics of such singular flows 
is the key in the analysis of general fluid flows governed by the well-known Euler 
(or Navier-Stokes) systems or their variants for both compressible and incompress- 
ible fluids, and has been one of the main focuses for applied analysists for decades. 
Substantial progress has been made in the past in studying singular flows by either 
rigorous analysis, or numerical simulations, or asymptotic methods [18]. In particu- 
lar, a rather complete theory exists for the 1-D shock wave problems, and both the 
theoretical understanding and numerical methods for 2-D smooth incompressible 
flows are quite satisfactory. Yet there are still many important issues to be settled, 
such as, the global existence of weak solutions and asymptotic structures of vari- 
ous approximate solutions generated by either physical considerations or numerical 
methods for such singular flows. 

One of the fundamental problems in the mathematical theory of shock waves 
for hyperbolic conservation laws is the well-posedness of the multi-dimensional gas- 
dynamical shock waves, for which the celebrated Glimm's method does not apply. 
Most of the previous studies along this line deal with either short time structural 
stability of basic fronts or asymptotic analysis and numerical simulations of dynam- 
ics of shock fronts. Due to the great complexity and the lack of understanding, it 
is reasonable for one to begin with some of the physically relevant wave patterns 
where a lot of experimental data, numerical simulations, and asymptotic results are 
available. One of the basic model for such studies is the supersonic flow past a 
pointed body [4], which is one of the fundamental problems in gas dynamics. We 
first study the stationary supersonic gas flow past an infinite curved and symmetric 
cone. The flow is governed by the potential equation as well as the boundary con- 
ditions on the shock and the surface of the body. This problem has been studied 
extensively by either physical experiments or numerical simulations. The rigorous 
analysis starts with the work of Courant and Friedrichs in [4] , where they show that 
if a supersonic flow hits a circular cone with axis being parallel to the velocity of 
the upstream flow and the vertex angle being less than a critical value, then there 
appears a circular conical shock attached at the tip of the cone, and the flow field 
between the shock front and the surface of the body can be determined by solving 
a boundary value problem of a system of ordinary differential equations. The local 
existence of supersonic flow past a pointed body has been established recently [1]. 
Our interest is on the structure of the global solution to this problem. We show 
that the solution to this problem exists globally in the whole space with a pointed 
curved shock attached at the tip of the cone and tends to a self-similar solution 
[2]. Our analysis is based on a global uniform weighted energy estimate for the 
linearized problem. The method we developed in [2] seems to be quite effective for 
other multi-dimensional problems. Indeed, similar approach can be used to study 
unsteady supersonic flows past a curved body for both potential flows and the full 
Euler system, and we obtain the local existence of shock waves in these cases [3] . 

Another challenging problem in the mathematical theory of fluid-dynamics 
is the theoretical foundation of the Prandtl's boundary layer theory [14]. In the 
presence of physical boundaries, the solutions to the invisicid Euler system cannot be 
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the uniform asymptotic ansatz of the corresponding Navicr-Stokes system for large 
Reynolds number due to the discrepancies between the no-slip boundary conditions 
for the Navier-Stokes system and the slip boundary condition for the Euler system. 
Indeed, the physical boundaries will creat vorticity and there is a thin layer (called 
boundary layer) in which the leading order approximation of the flow velocity is 
governed by the Prandtl's system for the boundary layers (see (3.1) in Section 
3). There are extensive literatures on the theoretical, numerical, and experimental 
aspects of the Prandtl's boundary layer theory [14]. Yet very little rigorous theory 
exists for the dynamical boundary layer behavior of the Navier-Stokes solutions for 
both compressible and incompressible fluids. One of the main difficulties is the 
well-posedness theory in some standard Holder or Sobolov spaces for the initial- 
boundary value problems for the Prandtl's system for boundary layers which is a 
severely degenerate parabolic-elliptic system, for which the only known existence 
results are proved locally in the analytic class in [15], except the series of important 
works of Olcinik who dealt with a class of monotonic data [12]. Indeed, Olcinik 
considered a plane unsteady flow of viscous incompressible fluid in the presence 
of an arbitrary injection and removal of the fluid across the boundary. Under the 
monotonicity assumption (see (3.5) in Section 3), Oleinik proved the well-posedness 
of local classical solutions to the Prandtl's system [12]. One of the open questions 
posed in [12] is to prove the global well-posedness of solution for the Prandtl's 
system under suitable conditions. Recently, we establish a global existence and 
uniqueness of weak solutions to the 2-D Prandtl's system for unsteady boundary 
layers in the class considered by Oleinik provided that the pressure is favorable. 
This is achieved by introducing a viscous splitting method and new weighted total 
variation estimate [16,17]. See Section 3 and [16,17] for more details. 

Finally, we turn to the motion of vortex sheets, which corresponds to a singular 
inviscid flow where the vorticity field is zero except on lower-dimensional surfaces, 
the sheets, and can be characterized as inviscid flows with finite local energy and 
with vorticity fields being finite Radon measures. The study on the existence and 
structure of solutions for the inviscid Euler system for incompressible fluids with 
data in such class is of fundamental importance both physically and mathemati- 
cally. Physically, vortex sheets can be used to model important flows such as high 
Reynolds number shear layers, and have many engineering applications. Mathemat- 
ically, the evolution of a vortex sheets gives a classical example of ill-posed problem 
in the sense of Hadamard, a curvature singularity developcs in finite time, and the 
nature of the solution past singularity formation is of great interest to know. This 
gives rise to many interesting yet challenging problems. Some of these are: Is there 
a well-posedness theory of classical weak solutions to the inviscid Euler system with 
general vortex sheets initial data? What are the structures of the approximate solu- 
tions to vortex sheets motions generated by either Navier-Stokes solutions or pratical 
numerical methods (such as particle method)? Can vorticity concentration and en- 
ergy defects occur dynamically? etc.. Despite the importance of these problems and 
past intensive effort in rigorous mathematical analysis, these problems are far from 
being solved. Better understanding has been achieved before singularity formation 
in the analytical setting, and studies on global weak solutions and their approxima- 
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tions start with the important works of Dipcrna-Majda [6]. Delort observes that no 
vorticity concentration implies that a weak limit in L 2 of an approximate solution 
sequences is in fact a classical weak solution to the two-dimensional Euler system, 
and thus proved the first existence of global (in time) classical weak solution to the 
2-D incompressible Euler equations with vortex sheets initial data provided that the 
initial vorticity is of distinguished sign [5]. Similar ideas have been used to study 
the convergence of approximate solutions generated by either viscous regularization 
[11] or partical methods [8,9] for vortex sheets with one sign vorticity. In the case 
that vorticity may change sign, the vortex sheets motion becomes extremely com- 
plex after singularity formation. Indeed, many important features of irregular flows 
seem to be connected with interactions and intertwining of regions of both positive 
and negative vorticities. Here we consider a mirror-symmetric flow which allows 
interactions but excludes intertwining of regions of distinguish vorticity, prove that 
there is no vorticity concentrations for approximate solution of such flow, and thus 
give the first global (in time) existence of vortex sheets motion with two-sign vor- 
ticity. We also consider the 3-D axisymmetric vortex motions. It is well-known 
that for smooth flows, the analysis for the axisymmetric 3-D Euler system without 
swirls is almost identitial to that of 2-D incompressible Euler system. However, 
this parallelness breaks down for the vortex sheets motions. Indeed, we will show 
that there exist no energy concentrations in approximate solutions to vortex sheets 
motion with one-sign vorticity for steady axisymmetric 3-D Euler system without 
swirls [7] . Some partial results on unsteady axisymmetric vortex sheets motion will 
also be discussed. 



2. The supersonic flow past a pointed body 

A projectile moving in the air with supersonic speed, is governed by the inviscid 
compressible Euler systems 

( d t p + div(pu) = 0, , 2 ^ 

where p,u — {u\, U2, U3) and p stand for the density, the velocity and the pressure 
respectively. We will only treat the polytropic gases so that p = p(p) = Ap 1 with 
gas constant A > and 1 < 7 < 3, 7 being the adiabatic exponent. 

Suppose that there is a uniform supersonic flow [u\, 112,113) — (0, 0,(70) with 
constant density po > which comes from negative infinity. Then the flow can 
be described by the steady Euler system. If we assume further that the flow is 
irrotational, so that one can introduce a potential function <f> such that u = V<f>. 

Then the Bernoulli's law implies that p = hr x ^q 2 + h(p ) - ±|V$| 2 ) = ff(V$), 

where h(p) is the specific enthalpy defined by h (p) = In this case, (2.1) is 

reduced to a second order quasilinear equation 

diw(iJ(V$)V$) = 0, (2.2) 

which can be verified to be strictly hyperbolic with respect to X3 if > c with c 
being the sound speed given by c 2 (p) = p (p). The flow hits a point body, whose 
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surface is denoted by m{x\, x 2 , £3) = 0. Since no flow can cross the boundary, the 
natural boundary condition is 

u ■ Vm = V$ • Vm = on m(xi, x 2l x 3 ) — 0. (2.3) 

If the vertex angle of the tangential cone of the pointed body is less than a 
critical value, it is then expected that a shock front is attached at the tip of the 
pointed body. Denote by fj,(xi, x 2 , £3) = the equation of the shock front, then the 
Rankine-Hugoniot conditions become 

• [i?(V$)V$] = 0, $ is continuous, on fj,(xi, X2, X3) = 0. (2.4) 

Our aim is to find a solution to this free boundary value problem, (2.2)-(2.4). 
When the pointed body is small perturbation of a circular cone, the local existence 
of solution to the problem (2.2)-(2.4) has been established in [1]. Our main goal is 
to establish a global solution. However, such a global shock wave might not exist 
in general for arbitrary pointed body due to the possibility of development of new 
shock waves in the large. Thus we assume further that the pointed body is a curved 
and symmetric cone. In this case, it will be more convenient to rewrite the problem 
(2.2)-(2.4) in terms of polar coordinates (r, 6, z) with r = \J x\ + x\ and z = X3. 
Assume that the tip of the pointed body locates at the origin, the equation of the 
surface of the body is r = b(z) with 6(0) = 0, and the equation of the shock front is 
r = S{z) with 5(0) = 0. Set $ = q z + ip(r, z). Then (2.2)- (2.4) become 

2 

((®, + d z ^f - c 2 ) d 2 zz 4>+ {{dr 4>) 2 - c 2 ) d 2 rr <p+2d r <p( qo +d z ^R 2 Z v - -d r <p = o, 

(2.5) 

-(go + dz <p)b (z) + d r tp = on r = b(z), (2.6) 

-[{q + d z ip)H]S'(z) + [d r ipH] = on r = S(z). (2.7) 

Moreover, the potential tp(r, z) is continuous on the shock, so it should satisfy 
(p(S(z),z) = 0. Then in [2], we have shown that the problem (2.5)-(2.7) has a 
globally defined solution as summarized in the following theorem: 

Theorem 2.1 Assume that a curved and symmetric cone is given such that 6(0) = 
0, 6'(0) = 6 , 6( fe )(0) = 0, 2 < k < h, and 

\z k -^(b(z)-b z)\<e Q for < k < k 2 , z > (2.8) 

with k\ and k 2 being some suitable integers. Suppose that a supersonic polytropic 
flow parallel to the z-axis comes from negative infinite with velocity qo, and density 
Po > 0. Then for suitably small Eq, 60 and g^ 1 , the boundary value problem (2.5)- 
( 2. 7) admits a global weak entropy solution with a pointed shock front attached at 
the origin. Moreover, the location of the shock front and the flow field between the 
shock and the surface of the body tend to the corresponding ones for the flow past 
the unperturbed circular cone r = b n z with the rate z~i. 
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It should be noted that there are no other discontinuities in our solution besides 
the main shock. Since the deviation of the surface of the body from that of a circular 
cone is sufficiently small (see (2.8)), any possible compression of the flow will be 
absorbed by the main shock. This is the mechanism to prevent the formation of 
any new shocks inside the flow field caused by the perturbation of the body. In 
particular, our results demonstrate that self-similar solution with a strong shock is 
structurally stable in a global sense. Indeed, the key element in the proof of Theorem 
2.1 is to establish some global uniform weighted energy estimates for the linearized 
problem of (2.5)-(2.7) around the self-similar solution with a strong shock obtained 
when the pointed body is a circular cone. This is achieved by a deliberate choice 
of multipliers which must satisfy a system of ordinary differential inequalities with 
complicated coefficients due to the structure of the background self-similar solution 
and the requirement of obtaining global estimates independent of z for the potential 
function and its derivatives on the boundary as well as its interior of a domain. 

When the projectile changes its speed, or it confronts some airstream, then 
the flow around the projectile will be time-dependent. Thus, we also consider the 
unsteady supersonic flow past a pointed body. Although our analysis applies to more 
general case [3], here we will only present the result for two-dimensional polytropic, 
unsteady, and irrotational flow past a curved wedge. For simplicity in presentation, 
we also assume that both the wedge and the perturbed incoming flow from infinity 
are symmetric about xi-axis. Let x 2 = b(xi) with 6(0) = be the equation of the 
wedge, and x 2 = S(t,xi) with S(t,xi = 0) = be the equation of the shock front. 
Then in terms of the velocity potential function <j) (so that u\ — d\4> and u 2 = d 2 4>), 
we are looking for solutions to the following initial boundary value problem 

dt (h{-4h + i|V0| 2 )) + Y? i=l d Xi (d Xi <f>H(-dt - i|V0| 2 )) = 0, 

t>0,x 2 >b 1 (x 1 ), (2.9) 
V<t>- (6'(ari),-l) = on x 2 =b(x 1 ), (2.10) 
[H]d t S+[Hd xl y\d xl S-[Hd X2 4>]=Q on x 2 =S(t, Xl ), (2.11) 

5(0, xi) = So(xi), (p(0,xi,x 2 ) = (p (xi,x 2 ), d t (j>(0,xi,x 2 ) = 0, (2.12) 

where Sq(xi) and 4>o(xi,x 2 ) are suitably small perturbations of the corresponding 
shock location and potential function respectively for the steady flow. Then we 
have the following local existence result [3] . 

Theorem 2.2 There exist positive constant Si and 5 2 and functions S(t,x\) and 
(j)(t, xi,x 2 ) defined on the regions {(t, x\)\0 < t < 5\, < x\ < 5 2 } and {(t, x\,x 2 )\0 
< t < 5\,0 < x\ < S 2 ,b(xi) < x 2 < S(t, x\)} respectively such that (S(t,x\), 
ip(t,Xi,x 2 )) solves the problem (2.9)-(2.12). 
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3. Prandtl's system for boundary layers 

Consider a plane unsteady flow of viscous incompressible fluid in the presence 
of an arbitrary injection and removal of the fluid across the boundaries. In this 
case, the corresponding Prandtl's system takes the form 

d t u + ud x u + v d y u + d x p = vdyyU, d x u + d y v = 0, (3.1) 

in the region, R = {(x,y,t)\0 < x < L, < y < +oo, < t < T}, where v, L and 
T are positive constants. The initial and boundary conditions can be imposed as 

( u\ t =o = u (x,y), u\ x = = ui(y,t), u\ y=Q = 0, v\ y = = v Q (x,t), . . 

\ \im y ^ +OQ u(x,y,t) = U(x,t), 

with U — U(x,t) given as is determined by the corresponding Euler flow. The 
pressure p = p(x,t) in (3.1) is determined by the Bernoulli's law: d t U + U d x U + 
d x p = 0. 

It follows from the physical ground that one may assume that 

U(x,t) > 0, u {x,y) > 0, ui(y,t) > 0, and v (x,t) < 0. (3.3) 

Due to the degeneracy in the Prandtl's system (3.1), the problem of well- 
posedness theory of solutions to the problem (3.1)-(3.3) in the standard Holder 
space or Sobolov space is quite difficult. In a series of important works by Olcinik 
and her coauthors [12], they studied this problem under the additional assumption 
that the data are monotonic in the sense that 

d y u Q {x,y) > 0, and d y Ui(y,t) > 0, (3.4) 

and prove that there exists a unique local classical smooth solution to the initial- 
boundary value problem (3.1)-(3.3) provided that the data are monotonic in the 
sense of (3.4). Here by local we mean that T is small if L is given and fixed, 
and T is arbitrary if L is small. One of the open problem in [12] is: What are 
the conditions ensuring the global in time existence and uniqueness of solutions to 
the problem (3.1)-(3.4) for arbitrarily given LI In [16,17], we study such problem 
and establish the global (in time) existence and uniqueness of weak solutions to the 
initial-boundary value problem (3.1)-(3.4) in the case that the pressure is favorable, 
i.e., 

d x p(x,t) < 0, t > 0, < x < L. (3.5) 
More precisely, we have ([16,17]): 

Theorem 3.1 Consider the initial-boundary value problem for the 2-D Prandtl's 
system, (3.1)-(3.2). Assume that the initial and boundary data satisfy the con- 
straints (3.3), (3.4) and (3.5). Then there exists a unique global bounded weak 
solutions to the initial-boundary value problem (3.1)-(3.2). Furthermore, these so- 
lutions are Lipschitz continuous in both space and time. 
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We remark here that the condition (3.5), that the pressure is favorable, is 
exactly what fluid-dynamists believe for the stability of a laminar boundary layer. 
This is also consistent with the case of stationary flows [12]. In the case of pressure 
adverse, i.e., ff > 0, separation of boundary layer may occur, so one would not 
expect the long time existence of solution to (3.1)-(3.3). Finally, we remark that in 
the case that (3.5) fails as in many pratical physical situations, short time existence 
of regular solution is still expected, which has not been established yet. 

4. Vortex sheets motions 

Two-dimensional vortex sheets motion corresponds to an inviscid flow whose 
vorticity is zero except on one-dimensional curves, the sheets. Thus, it is governed 
by the following Cauchy problem 

J dt w + u ■ Vw = 0, u = K * w, . . 

\ w(t = 7 x)=w (x)eM{R 2 )nH lo 1 c (R 2 ), (4 ' lj 

where u — (ui,u 2 ) is the velocity field, w = V 1 - ■ u is the vorticity, and K is the 
Biot-Sawart kernel, and Ai(M. 2 ) denotes the space of finite Radon measures in M 2 . 
Let (u e ,w e )(w £ = V -1 • u e ) be a sequence of approximate solution to (4.1) with the 
following bounds 

sup / \w 2 (x,t)\dx < Ci(T), sup / \u £ (x,t)\ 2 dx < C 2 {T,R). (4.2) 

0<t<T JR2 0<t<T J\x\<R 

Then there exist ueL°° (0,T 7 L 2 oc (M. 2 )) and weL°° (0,T 7 M(R 2 )) with w = 
V 1 - ■ u such that 

w £ ^winM ([0, T] xl 2 ), vl e — ^ u in L 2 loc ([0, T] x E 2 ) (4.3) 

Two main questions arise: Is (u, w) a classical weak solution to problem (4.1)? 
Does either vorticity concentration (i.e. lim e r ^ + Jb \w £ (x,t)\dx > 0) or energy 
defects ( lim e ^ J\ x \<r \ u£ ( x , t)\ 2 dx > Ji x i <r \u(x, t)\ 2 dx for some t > 0, R > 
0) occur? It is clear from the structure of the 2-D Euler system that no energy 
defects implies strong L 2 -convergence of the velocity field and thus the existence 
of the classical weak solution to (4.1). A less obvious fact is that no vorticity 
concentration also implies the weak limit (u,w) being a classical solution to (4.1), 
which follows from the vorticity formulation of (4.1) as observed by Delort [5]. 
Then the convergence to a classical weak solution to the Cauchy problem (4.1) is 
proved for approximate solutions generated by either regularizing the initial data 
[5], or Navier-Stokes approximation [11], or vortex blob methods [8], or point vortex 
methods [9] provided that the initial vorticity is of distinguished sign. To study the 
corresponding issues for flows where interactions of regions of both positive and 
negative vorticities are allowed, in [10], we study a 2-D mirror-symmetric vortex 
sheets motion whose vorticity is a integrable perturbation of a non-negative mirror- 
symmetric radon measure. Here a Radon measure fi is said to be non-negative 
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mirror-symmetric (NMS) if ^| R 2_ > and \i is odd with respect to x\ = 0. Then 
the main results in [10] can be summarized as 

Theorem 4.1 Assume that w = m + [i 2 such that nieM c (R 2 ) n if z ~*(R 2 ) n 
NMS and /12 e L;i(R 2 ). Then there exists a global (in time) classical solution 
(u,w)eL°° (0,T,L 2 (R 2 ))®L°° (0, T,M(R 2 )) to the initial value problem (4.1). Fur- 
thermore, this weak solution can be obtained as a limit of either a sequence of smooth 
inviscid solutions or a sequence of solutions to the Navier-Stokes system. 

It should be noted that this is the only result of existence of classical weak 
solution to (4.1) involving vorticities with different signs. This is proved by showing 
Jo su -Px eRi S B {x ,6) \w £ {y,t)\dydt -> as 6 -> 0+ uniformly in e, i.e., no vor- 
ticity concentration occur any where. Theorem 4.1 also indicates that interactions 
of regions of positive and negative vorticities without interwining may not cause 
concentration in vorticity. There remain many important open problems for the 
2-D vortex sheets motion such as the existence of classical weak solution to (4.1) 
for general vortex sheets initial data, and whether energy defects occur dynamically 
even in the case of one-sign vorticity. 

Finally, we consider 3-D axisymmetric vortex sheets motions. In cylindrical 
coordinate, (r, 8, z), axisymmetric solutions of 3-D Euler system have the form 

u(x, t) = u r (r, z, t)e r + u e (r, z, t) eg + u z (r, z, t) e z , p(x, t) = p(r, z, t) (4.4) 

where e r = (cos 8, sin 6, 0), eg = (— sin#, cosd, 0), and e z = (0,0,1). The axisym- 
metric flow is said without swirls if u e = 0. In this case, the vorticity field is 
given by w = V x u = w e e e with w e — d r u z — d z u r , and _D t (r _1 w e ) — with 
Dt = d t + u r d r + u z d z . Thus, for smooth data, the theory for 3-D axisymmetric Eu- 
ler system without swirls is almost parallel to that of 2-D Euler equation. However, 
this similarity breaks down for vortex sheets motions. Indeed, in sharp contrast to 
the 2-D case, we show in [7] that there are no energy defects for a suitable sequence 
of approximate solutions for 3-D steady axisymmetric vortex sheets motion with 
one-signed vorticity. Precisely, we have 

Theorem 4.2 Let (u e ,p e ) be smooth axisymmetric solution to the 3-D steady Euler 
system: {u E ■ V)u £ + Vp e = f £ , divu £ = 0, ieK 3 , for some given axisymmetric 
function f e with f e — / weakly in L 1 (M 3 ). Suppose further that 



where w £ — V x u £ = (w £ ) eg. Let u be the weak limit of u £ in L 2 (R 3 ). Then u is 
a classical weak solution to (u ■ V)m + Vp = /, divu = 0, xeR 3 . Moreover, there 
exists a subsequence {u £j } of {u £ } such that u £j converges to u strongly in L 2 (R 3 ). 

The proof of this theorem is based on a shielding method and the following 
fact which is valid for both steady and unsteady flows [7]. 




(4.5) 
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Theorem 4.3 Let {u £ } be a sequence of approximate solutions for 3-D axisymmet- 
ric Euler system with general vortex-sheets data generated by either smoothing the 
initial data or Navier-Stokes approximations. Let A be the support of the defect mea- 
sure associated with u e eL 2 ([0,T] x M 3 ). If A ^ <p, then A n{(x,t)eM. 3 x [0,T]| r > 

Theorem 4.3 implies that if there arc no energy defects away from the sym- 
metry axis, then strong L 2 -convergence takes place. It remains to study whether 
energy defects occur for unsteady axisymmetric 3-D Euler equations. 
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